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Experimentally,  the  investigation  of  reactions  in  explfe- 
sive  materials  has,  until  fairly  recently,  been  concerned  largely  with 
fast  reactions.  There  are  obvious  difficulties  inherent  in  these 
experiments  due  to  the  very  short  time  in  which  processes  occur,  the 
complication  of  mass  flow  (that  is  massive  transport  of  material  as, 
for  example,  in  the  region  of  a  shock  or  pressure  wave)  and  the 
inability  to  determine  the  thermodynamic  parameters  due  to  extreme 
temperatures  and  pressures.  These  experiments  have,  however,  led  to 
good  phenomenological  descriptions  of  the  propagation  of  detonation. 
However,  -they  have  yielded  little  insight  into  the  mechanism  of  growth 
of  reaction,  particularly  in  regard  to  the  nature  of  reaction  at  the 
molecular  level.  A  knowledge  of  the  conditions  necessary  for 
initiation  of  detonation  is  of  great  practical  importance  both  from 
the  viewpoint  of  causing  a  detonation  to  occur  when  it  is  desired  and 
also  to  prevent  it  from  occurring  when  it’ is  not.  This  problem  has 
received  considerable  study  and  many  methods  have  been  designed  for 
its  achievement.  However,  since  the  detailed  mechanism  of  initiation 
at  the  molecular  level  is  unknown,  this  exact  control  of  initiation 
is  not  presently. attainable. 


It  seems  likely  that  the  stimulus  responsible  for  the 
subsequent  release  of  chemical  energy  is  primarily  thermal  and  that 
other  practical  means  of  producing  initiation  (the  application  of 
shock,  mechanical  action,  electrostatic  charge,  sparks  and  high  inten¬ 
sity  light  pulses)  are  probably  thermal  in  their  effect.  An 
investigation  was  therefore  undertaken  of  the  thermal  reactions  in 
solids  which  may  subsequently  grow  into  deflagration  or  detonation. 

•  Scientifically  initiation  is  an  indeterminate  property,  however,  it 
includes  the  point  in  time  at  which  the  growth  reactions  become 
apparent;  Thepe  growth  reactions  can  be  classified  into  three 
general  kinds: 

a.  The  first  involves  "chains"  of  reactions  within  the 
-  bulk  of  the  *solid  or  at  the  surface  of'  the  solid,  the  chain  being 
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propagated  by  either  products  of  the  reactions  (atoms  in  new  charge 
states,  charge  carriers  themselves,  or  electronic  excitation)  or 
vibrational  excitation,  where  the  energy  released  by  individual  steps 
in  the  process  may  be  considered  to  be  contributing  to  local  changes 
in  temperatures  ("hot  spot"  model). 

b.  The  second  type  deals  with  processes  within  the  material 
leading  to  products  which  are  transported  to  sites  on  the  surface  or 
at  grain  boundaries  and  the  like.  The  products  would  be  then  consid¬ 
ered  to  aggregate  to  form  new  unstable  phases. 

c.  The  third  category  considers  the  possibility  of  a  direct 
phase  change  of  the  material  accompanied  by  a  release  of  energy,  the 
material  having  become  vibrationally  unstable.  This  latter  type  of 
model  differs  importantly  from  the  first  two  classes  in  that  no  trans¬ 
port  of  reaction  products  is  involved.  Although  this  mechanism  does 
not  account  for  the  known  effect  of  relatively  small  quantities  of 
impurities  and  the  apparent  importance  of  surface  reactions  in  explo¬ 
sives,  a  relatively  large  number  of  phase  changes  have  been  observed 
in  explosive  materials  and  the  possibility  of  a  large  energy  release 
after  a  phase  change  is  of  interest. 

In  the  study  of  these  various  models  of  growth  reactions, 
the  technique  of  neutron  scattering  can  provide  information  which  is 
essential-  to  the  understanding  of  the  processes  taking  place  at  the 
molecular  level.  These  processes  involve  both  intra-  and  inter-molec¬ 
ular  vibrations  of  the  molecular  arrangement  and  their  changes  at 
different  stages  of  the  reaction.  The  third  mechanism  could  be  in¬ 
vestigated  by  noting,  for  example,  that  lattice  vibrational  modes 
must  show  the  onset  of  a  phase  change.  The  results  of  thermal 
decomposition  or  irradiation  experiments  or  the  known  effects  of 
imperfections  (dislocations,  defects  or  impurities)  on  the  explosive 
sensitivity  can  now  be  investigated  on  the  microscopic  level.  An 
attempt  to  correlate  these  findings  with  the  microscopic  properties 
of  explosive  or  related  materials  may  therefore  become  possible. 

This  paper  reports  the  preliminary  results  of  these  studies. 
First, 'we  Shall  describe "the  ways  in  which  .neutrons  may-  interact  with 
nuclei  in  a  solid  or  liquid  and  the  various  types  of  information 
which  can  be  obtained  from  these  measurements . 

Neutron  diffraction  is  now  used  quite  extensively,  often  as 
a  necessary  complement  of  x-rays,  to  determine  the  position  of  light 
atoms  (especially  hydrogen)  in  a  crystalline  structure.  Neutrons  of 
thermal  energy  (kT=25.  millielectron  volts)  are  scattered  elastically 
(without  exchanging  energy)  and  the  neutron  waves  scattered  from  the 
periodically  arranged  nuclei  give  rise  to  interference  effects  or 
Bragg  reflection.  We  obtain,  in  this  case,  a  "static"  picture  of  the 
crystal,  that  is,  every  atom  is  vibrating  around  its  equilibrium. 

•  position.  Neutrons  can  also  be  used  as  "probes"  to  study  the  dynamics 
of  motions  of  atoms  and  molecules  in  the  crystal.  If  their  incident 
energy  is  small  (5  mev)  compared  to  the  vibrational  energy  of  the 
atom  ip  the  crystal  (kT=25  mev),  energy  will  be  gained  by  these 
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neutrons  in  a  scattering  event  and  an  analysis  of  the  distribution  in‘ 
energy  of  these  neutrons  after  scattering  will  reflect  the  dynamics 
of  motion  inside  the  crystal.  Those  motions  are  observed  which  take 
place  in  the  far  infrared  frequency  region  (from  40  to  800  cm”1) 
without  being  subject  to  symmetry  dependant  selection  rules,  so  that 
all  possible  transitions  are  observable.  Although  the  energy  resol¬ 
ution  is  not,  as  yet,  comparable  to  optical  spectroscopic  methods, 
the  technique  of  neutron  inelastic  scattering  offers  a  complement  to 
Infrared  or  Raman  techniques .  It  should  also  be  noted  that  this 
technique  is  most  sensitive  to  the  motions  of  hydrogen  atoms  since 
the  scattering  cross-section  for  hydrogen  is  considerably  larger 
80  barns)  than  for  any  other  atom.  In  the  case  of  monatomic 
solids  with  cubic  symmetry,  the  frequency  distribution,  g(y),  (i*e., 
the  number  of  states  occurring  per  unit  frequency  interval)  can  be 
obtained  exactly  from  the  measured  neutron  differential  inelastic 
scattering  cross-section.  Only  approximate  values  of  g(^))  can  be 
obtained  in  the  case  of  polyatomic  molecules  with  a  symmetry  other 
than  cubic.  From  the  frequency  distribution  fundamental  macroscopic 
properties  such  as  specific  heat  (Cv),  entropy  (S),  internal  energy 
(E)  and  free  energy  (F)  can  be  derived. 

In  this  work,,  vibrational  frequencies  in  the  far  infrared 
region  (30-1000  cm”1;)  which  have  not  been  obtained  previously  or  are 
optically  inactive  will  be  presented  for^  and  Y  HMX,  RDX,  picric 
acid  and  lead  styphnate.  This  frequency  region  is  of  importance 
since  it  is  expected  to  contain  most  lattice  modes  (inter-molecular 
vibrations)  and  also  because  the  low  intra-molecular  frequencies  are 
more  sensitive  than  the  higher  frequency  vibrations  to  impurities  in 
the  lattice  or  to  change  in  molecular  conformation.  The  derived 
frequency  distribution  wili  be  used  to  calculate  thermodynamic  quan¬ 
tities  and  to  compare  them  with  measured  values  where  available. 

$  and  HMX  and  RDX:  These  nitramine  compounds  are 
secondary  explosives  with  comparable  sensitivity.  RDX  forms  a  six 
membered  ring,  whereas  HMX  forms  an  eight  membered  ring  and  can  exist 
in  several  polymorphic  forms  > /£ ,  Y  ,  J  )  -  Of  the  two  forms 
which  have  been  studied  in  this  work,  /£  HMX  forms  the  more  compact 
structure,  the  ^  form  being  more  open .  < 

The  spectra  of  scattered  neutrons  for  these  three  explo¬ 
sives  are  shown  in  Figures  -1  and  2 .  Table  1  gives  -a  summary  of  the 
inelastic  neutron  scattering  results  for  RDX  and  HMX  and  tentative 
assignments  of  the  vibrational  frequencies  together  with  available 
infrared  data  of  these  compounds  in  the  frequency  range  below  1200 
cm'1.'  The  similarity  of  the  vibrational  frequencies  in  these  com¬ 
pounds  is  not  unexpected  since  they  have  the  same  chemical  groups 
and  comparable  bond  distances  and  bond  angles.  It  should  be  noted 
that  in  the  case  of  ft  and  yUMX,  the  neutron  spectra  show  differences 
in- shape,  especially  in  the  region  of  lattice  modes  (an  additional 
peak  at  105  -cm"1  is  observed  in  ^HMX).  These  spectral  differences 
show  the  sensitivity  of  the  neutron  scattering  technique  to  confor¬ 
mational  changes  and  differences  in  intermolecular  crystalline  forces. 
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The  frequency  distribution  obtained  from  the  neutron  spectra  of  com¬ 
plex  molecules  has  the  greatest  validity  .in  the  region  of  lattice 
modes,  since  for  these  modes  the  molecules  move  as  single  entities 
similar  to  atoms  in  a  monatomic  crystal.. 

The  calculation  of  Cy  vs  T  has  been  performed  in  the  case 
of  RDX  and  ft  and  ^HMX,  using  the  frequency  distributions  g(\^), 
up  to  500  cnT^,  derived  from  the  time-of-f light  data.  The  low  tem¬ 
perature  portion  of  this  curve  offers  the  most  reliable  values  of  Cy 
since  lattice  modes  represent  the  major  contribution  to  Cv  in  that 
region.  It  is  unfortunate  that  in  the  case  of  HMX  there  are  no 
measured  values  of  the  specific  heat  to  compare  with  these  results. 
For  RDX,  Cp  has  been  measured  only  in  a  narrow  temperature  range 
(around  room  temperature) .  The  change  in  free  energy  A  F,  internal 
energy -A  E  and  entropy  A  S  for  the  transf ormation  n  HMX  HMX 

are  estimated  to  be  +1  Kcal/mole,  +0.144  Kcal/mole  and  +2.8  entropy 
units,  respectively.  The  A  E  qf  this  transformation  is  very  small, 
but  this  is  not  unexpected  since  no  bond  is  broken  in  the  transfor¬ 
mation.  Only  a  change  in  the  intermolecular  forces  takes  place,  due 
to  different  packing  and  molecular  environment.  The  inprease  in 
entropy  is  going  from  /to  Y  HMX  suggests  that  the  J  structure  is 
more  disordered  than  the  structure.  Neutron  diffraction  data  of 
HMX  obtained  in  this  laboratory  have  shown  a  certain  amount  of 
rotational  freedom  of  the  CH2  groups  in  the  lattice.  These  CH2 
librational  motions  (hindered  rotations)  probably  increase  in  ampli¬ 
tude  in  the  /  s truer urn. 

Lead  Styicr-ite,  Picric  Acid,  Stvphnic  Acid:  The  study 
of  these  three  compounds  was  undertaken  because  they  provide  a  series 
.having  similar  molecular  structure  but  vastly  different  explosive 
properties.  Lead  ssyphnate  is  one  of  the  most  sensitive  primary 
explosives  known,  ccvpatable  to  lead  azide  and  nitroglycerin  in 
impact  sensitivity.  Picric  acid  is  also  an  explosive  but  much  less 
sensitive  to  impact  than  lead  sbyphnate,  while  styphnic  acid  is  the. 
most  stable.  There  ha.i  n<-en  considerable  discussion  concerning  the 
still  unresolved  structure  of  normal  lead  sty.phnate.  Two  structures 
have  been  proposed: 
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,  Because  no  large  single  crystals  of  normal  lead  styphnate 

have  thus  far,  been  available,  it  has  not  been  possible  to  employ 
x-ray  or  neutron  diffraction  techniques  to  determine  the  structure 
directly.  Infrared  measurements  have  been  performed  in  the  600-4000 
cm"1  frequency  range\?hich,  because  of  the  absence  of  a  peak  in  the 
3100-3400  ctn’*  region,  have  been  interpreted  to  suggest  the  absence 
of  water  of  hydration.  However,  no  line  corresponding  to  0-H 
stretching  is  observed  either  so  that  the  results  are  inconclusive. 

A  NMR  study  of  lead  styphnate  in  H2O  solution  also  failed  to  identify 
the  correct  structure? 

It  was  hoped  that  a  study  of  the  neutron  spectrum  of  hydra¬ 
ted  lead  styphnate  could  give  evidence  as  to  whether  I  or  II  was  the 
correct  structure.  The  neutron  spectrum  of  lead  styphnate  is  shown 
in  Figure  3.  This  spectrum  is  not  too  different  from  those  obtained 
from  a  large  number  of  hydrates.  The  intense  peaks  at  430  and  350 
cm"1  might  correspond  to  a  composite  of  transitions  due  to  H2O 
■  librations  and  Pb-0H2  stretching  and  those  at  230  and  200  cm**1  might 
be  due  to  hydrogen  bond  stretching-,  all  of  which  are  observed  in 
other  hydrates.  Unfortunately,  the  interpretation  of  this  spectrum 
is  complicated  by  the  preseitee  of  low  energy  intramolecular  modes  of 
the  styphnate  ring. 

.  The  spectra  of  neutrons  scattered  from  styphnic  acid  and 
picric  acid  are  shown  in  Figure  A (a)  and  4(b).  Because  of  the 
similarity  of  the  structure  of  chese  three  molecules,  one  might  be 
able  to  identify  the  intramolecular  vibrations  of  the  ring  itself  if 
they  are  common  to  all.  Table  2  summarizes  the  neutron  scattering 
data  for  these  three  compounds,  tentative  assignments  of  the  vibra¬ 
tional  frequencies  and  also  the  available  infrared  data  for  the 
•  frequency- range  below  1200  cm”1-.  In  styphnic  acid  in  which  no  water 
of  hydration  exists,  a  peak  is  observed  at  415  cm"1,  while  none  occurs 
in  picric  acid.  This  mode  of  vibration  may  also  give  rise  to  the 
430  cm”1  peak  in  lead  styphnate.  Consequently,  although  the  neutron 
spectrum  of  normal  .lead  styphnate  is  very  similar  to  spectra  of  salts 
containing  water  of  hydration,  ‘the  existence  of  ring  deformation 
’  frequencies  at  430  and  350  cm"1  in  both ‘the  neutron  spectra  of  lead 
styphnate  and  styphnic  acid  seem  to  favor  structure  II  of  lead 
styphnate. 


The  general  features  of  the  spectra  of  these  three  compounds 
which  have  greatly  different  physical  properties  are  of  considerable  ' 
interest.  The  neutron  spectrum  of  picric  acid,  a  mild  explosive,  is 
quite  complex,  displaying  several  sharp  peaks  in  the  frequency  range 
above  300  cm"1  and  these  are  probably  due  to  intramolecular  modes. 

The  spectrum  of  styphnic  acid,  which  is  not  an  explosive’,  shows  a 
continuous  increase  in  scattered  neutron  intensity  from  the  region 
of  lattice  modes  ( A-'lOO  cm"1)  to  that  of  intramolecular  vibrations 
(  >  200  cm-1)  with  no  clear  separation.  The  spectrum  of  normal 
l'ead  styphnate  which  is  the  most  sensitive  to  impact  shows  some 
separation  between  inter-  and  intra-molecular  modes  with  intensity 
considerably  lower  in  the  region  of  intermolecular  modes  or 
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lattice  vibrations.  To  quantitatively  interpret  these  differences  in- 
terms  of  differences  in  inter-  and  intra-molecular  forces,  and  to 
relate  them  to  thermodynamic  and  explosive  properties  would  require 
a)  a  calculation  of  the  vibrational  frequencies  of  the  isolated  mole¬ 
cules  and  this  has  been  undertaken  and  b)  a  detailed  knowledge  of  the 
crystal  structures.  However,  good  single  crystals  of  these  materials 
have  not,  as  yet,  been  obtained.  Finally,  the  values  of  Cv  as  a 
function  of  temperature  for  picric  acid  and  lead  styphnate  have  been 
calculated  from  the  neutron  data .  Only  values  of  Cp  have  been 
measured  and  these  in  a  narrow  temperature  range  in  the  vicinity  of 
300°K. 


Correlation  with  Sensitivity;  In  addition  to  the  assign¬ 
ment  of  vibrational  frequencies  in  the  0-800  cm“l  frequency  range  an 
attempt  has  been  made  to  correlate  our  microscopic  observations  with 
some  of  the  physical  properties  of  the  explosives  under  study. 

In  a  molecular  crystal  there  are  not  only  intramolecular 
vibrations  (within  a  given  molecule)  but  also  intermolecular  motions 
in  which  adjacent  molecules  vibrate  relative  to  each  other.  These 
intermolecular  motions  (lattice  modes)  are  generally  of  lower  energy 
than  the  intramolecular  motions  since  the  forces  binding  adjacent 
molecules  are  weaker  than  those  binding  the  atoms  within  the  molecule. 
The  probability  of  thermally  exciting  a  vibration  whose  energy  is  £, , 
is  given  by  . 

P  ( c. )  =  (exp  -l)-l . 

An  examination  of  this  expression  shows  that  the  thermal  energy 
available  to  a  crystal  will  be  primarily  deposited  in  the  low  fre¬ 
quency  vibrations  which  the  crystal  can  sustain.  The  amount  of 
energy  given  to  the  intramolecular  vibrations  will  therefore  depend 
very  strongly  upon  the  relative  density  of  the  low  energy  lattice 
vibrations.  This  consideration  should  be  of  importance  since  the 
molecular  vibrations  strain  those  intramolecular  bonds  whose  rupture 
leads  to  a  release  of  large  amounts  of  energy.  In  general,  therefore, 
one  would  expect  that  the  larger  the  fraction  of  available  thermal 
energy  which  is  deposited  in' intramolecular  vibrations,  the  easier 
will  it  be  to  rupture  molecular  bonds. 

The  lattice  modes  also  influence  the  first  mechanism  of 
the  growth  reaction  described  earlier,  since  the  probability  of  for¬ 
mation-  of  a  "hot  spot"  in- a  localized  region  of  the  crystal  is- 
dependant  upon  the  ease  with  which  the  crystal  lattice  is  able  to 
dissipate  thermal  energy.  This  conductivity  is  determined  by  the 
lattice  vibrations.  Thus  lattice  modes  condition  the  crystal  to 
promote  the  formation  of  "hot  spots"  which  will  subsequently  lead  to 
explosion. 

An  attempt  has  therefore  been  made  to  relate  the  vibration¬ 
al  and  physical  properties  of  the  compounds  under  study.  The  inter¬ 
nal  energy  E  (T)  has  been  calculated  from  the  relation 

-  re2 

E(T)  -  3Nr  j£!  |  coth  •£  g(£ )  d  C 
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where  3Nr  is  the  number  of  degrees  of  freedom,  is  the  vibrational 
energy  and  g(£)  is  the  frequency  distribution  derived  from  the  neu¬ 
tron  data.  This  expression  is  rigorously  valid  in  the  harmonic 
approximation  only  and,  as  pointed  out  previously,  the  applicability 
of  gCo)  is  limited  by  various  approximations.  The  ratio  R  is  formed 
of  the  internal  energy  of  the  low  energy  intramolecular  vibrations 
(-200-800  cm"1)  to  the  internal  energy  of  the  lattice  modes  (0-200  cm"  i 
The  value  of  this  ratio  has  been  plotted  in  the  figure  shown  below 
against  the  measured  impact  sensitivity  for  each  of  the  four  explo¬ 
sives  considered  in  this  work.  The  impact  test  value  is  a  minimum 
height  at  which  at  least  one  of  ten  trials  results  in  explosion.  This 
value  for  each  material  depends  upon  the  test  conditions.  However, 
there  is  general  agreement  concerning  the  order  in  which  explosives 
should  be  classified.  ....  .  . 


Plot  a  ‘  Plot:  b 


Plot.ji  substantiates  the  qualitative  arguments  presented  above  and 
shows~clearly  that  ttie  explosive  becomes  more  sensitive  as  R  increases. 
The  ratio  R  was  also  plotted  as  a  function  of  temperature  (Plot  _b)  and 
according  to  this  calculation,  the  sensitivity  of  the  four  explosives 
should  vary  at  different  rates  with  temperature.  In  particular 
,/^HMX  and  picric  acid  should  show  the  greatest  difference  in  the 
rate  of  change  near  room  temperature  (/0HMX  having  the  slowest  rate 
and  picric  acid  the  greatest) .  Measurements  of  the  temperature 
dependance  of  the  impact  sensitivity  of  these  explosives  would 
constitute  a  very  good  test  of  these  calculations. 
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Finally  in  the  case  of  explosives  where  the  impact  sensi¬ 
tivity  is  known  to  vary  with  the  amount 'of  '“doping",  a  study  of  the 
corresponding  changes  in  the  value  of  the  ratio  R  with  the  amount  of 
impurities  will  be  undertaken  in  order  to  pro  ,’ide  additional  evidence 
substantiating  the  correlation  described  in  this  paper. 


\ 


CONCLUSION 


The  inelastic  scattering  of  neutrons  provides  a  means  of 
studying  explosives  by  which  the  nature  of  molecular  motions  can  be 
determined  and  through  which  molecular  dynamics  can  be  correlated 
with  macroscopic  properties  such  as  thermodynamic  quantities.  A 
tentative  correlation  has~beenmade  between  explosive  sensitivity 
and  internal  energy  for  four  explosives.  These-  results  show  the 
importance  of  understanding  the  microscopic  behavior  in  order  to  pre* 
diet  explosive  characteristics . I  - 


A  number  of  additional  experiments  are  planned  for  a  more 
complete  understanding  of  the  mechanism  of  the  growth  of  reactions 
at  the  molecular  level  and  to  further  test  the  correlation  obtained 
in  the  present  work. 
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Table  I 


Frequencies  of  Vibration 

(0-1200  cm“l) 

/Ihmx 

RDX  Assignments 

Neutron*- 

I.R.2 

Neutron*- 

I.R.2 

Neutron*-  I.R."* 

-  - 

1202 

1280 

1210 

1148 

. 

1090 

1088 

1035 

1015 

1015 

968 

. 

967 

949 

940 

919 

b-870 

.911 

b-890 

864 

862 

- - —  -  -  . 

.  832 

840 

• 

773 

764 

772 

762 

732 

751 

b4  -670 

707 

737 

b-670 

600  jring  defor- 

mat ion  or 


416 

424 

424 

CH  bending 

-  336 

330 

338  _ 

180 

'  235 

200  . 

NO2  torsicns 

120 

117  — 

'  104 

• 

80 

69 

70 

lattice 

■  .  .  • 

55 

modes 

40 

48 

35  ' _ 

1)  This  work. 

■ 

■ 

2)  -Bedard  et  al . , 

Can.  J.  Chem.  40, 

2278  (1962). 

3)  Hyers  &  Bedard,  CARDE  Tech  Memo  139/57,  Jan.  1957 


4)  Prefix  b  denotes  a  broad  peak. 
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Table  II 

Frequencies  of  Vibration  (0-1200  cm"^) 


Pb  Styphnate _ Styphnic  Acid _ Picric  Acid 


Neutron3 


430 

350 

230 


I.R.1 2 

Neutron3 4 

I.R.3 

Neutron3 

I.R. 

•  1180 

1130 

1100 

1080 

1080 

1080 

• 

980 

930 

940 

940 

-  -  - 

915 

860 

772 

775 

755 

735 

730 

775 

720 

718 

705 

690 

685 

690 

685 

— 

525 

« 

415 

330 

360 

230 

•t 

t; 


l 


) 

1 


200 

• 

193 

156 

148 

97 

b-105 

W 

74 

80 

64 

40 

48' 

48 

36 

30 

38 

28 

broad 

band 


1)  This  work. 

2)  R.  Zingaro,  J.  Amer.  Chem.  Soc.  76,  816  (1954). 

3)  F.  Pistera  et  al..  Anal.  Chem.  22,  495  (1960). 

4)  Prefix  b  denotes  a  broad  peak. 
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